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Abstract 



Nonperturbative corrections up to m^ 2 to the heavy lepton energy distri- 
butions are investigated in the inclusive semileptonic weak decays of heavy 
flavors in QCD. In the case of S-meson decays, for b — * utv transitions they 
decrease the decay rate by 6% of its perturbative value, while for b — > ctu 
they decrease it by 10%. 
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I. INTRODUCTION 



In this paper we investigate the power corrections to the heavy lepton energy distribution 
in the inclusive decays 

H b -> ruX, 

where Hb is a hadron containing heavy quark b and r is the r-lepton. We include non- 
perturbative corrections up to the order 1/m 2 , and show that they cause the decay rate of 
i?-mesons to decrease by 6 to 10 percent of its perturbative value depending on the mass of 
the quark in the final state. The general approach of this work is based on the heavy quark 
mass and the operator product expansions 0- @. The massless lepton case was analyzed 
in the works ||- [ I2"f . In those works the derivation of the corresponding matrix element 



and the operator product expansion were discussed in details. The lepton mass does not 
effect neither hadronic nor weak leptonic tensors. Therefore we can use expressions (l)-(6) 



and (Al)-(All) of Ref. [ I I for the hadronic invariant functions to get the corresponding 
matrix element. What differs the decay with a heavy lepton from the decay with a massless 
lepton is the phase space of the particles in the final state. The present paper is a continu- 



ation of our previous work fllT|. Here we will not repeat the discussion of applicability and 



interpretation of the results. Everything said about it in [JTTJ] is valid in the cases considered 
in this paper. 

The first experimental observations of the decay were made recently using the missing 
energy tag H, [|]. The missing energy was associated with two v T in the decay chain 
b — > t~v t X, r~ — > u T X', which made it difficult to reconstruct. The branching ratio 
was found to be 4.08 ± 0.76 ± 0.62% in and 2.76 ± 0.47 ± 0.43% in @ (more recent 
analysis), which is compatible with the Standard Model. However because of the difficulties 
in identification of the decay mode Hb — > tvX , the accuracy of the measurements is still 
insufficient for direct comparison of the results of this paper with the data. 



II. THE HEAVY LEPTON ENERGY DISTRIBUTIONS 



As it was discussed in [[LT| in the semileptonic inclusive decays there are three independent 
kinematical variables: E T , go, and g 2 , where E T is the energy of the emitted lepton (r in this 
paper), go is the energy and q 2 is the invariant mass of the lepton pair. The fully differential 
distribution can be written as follows: 

dT . 2m2 



dE T dq dq 2 128 tt 2 



\M(E T ,q ,dq 2 )\ 2 . (1) 



The matrix element of the process A4(E T ,q ,dq 2 ) involves the same hadronic and leptonic 
tensors as used in ]nj] and which are given by the expressions (l)-(6) and (Al)-(All) of 
Ref. The corresponding phase space and the boundaries for the kinematical variables 



for the massive lepton case are briefly discussed in Appendix. 

The double differential distribution in E T and q 2 could be obtained from ([!]) by inte- 
grating over go- The corresponding expressions are somewhat cumbersome because of the 
complicated limits of integration over g in (g 2 , g^) plane. Also there is little use of the 
double distribution because of the difficulties in observing the considered decays at present 
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time. Therefore we will not write out the formulae for the double distributions in (q 2 ,E T ), 
but will rather integrate Eq. ([!]) twice and get the energy distribution for the charged lepton 
in the final state. One more remark is in place here. Although we are looking at a hadron 
decay, in our dynamical approach we are considering a heavy quark decaying in the external 
field created by its interactions with the light degrees of freedom. That is why in the case 
at hand we have to use the quark kinematical boundaries rather then the hadronic ones. 
This means that in the formulas (|A7[)-(X9") we have to use mj and m q instead of Mn b and 
Md correspondingly. We will use such a kinematical boundary in the following calculations. 
Our results then should be understood in the sense of duality (see discussion in JTT[] ) . 

To get the energy distribution we integrate (|1|) over dq 2 and dq . In the dimensionless 
variables: 

2E T m\ ml 

x = , p q = — | , and p T = — 2 , (2) 

m b m£ mi 

the result of the integration then takes the following form: 

1 dT 



2^x 2 - 4p T {3(1 - p q + p T )x-2x 2 -Ap T + Pq (3-p q + 3 p T ) f+ 

(-3 - Pq + 6p T - p qPr - 3p 2 T ) f + 2 (1 - p T f f + 

5x 2 14 p 2 2 2 

K b [ — + ~ 3 Pi ( 3 " P, + 3pr) / + (1 + 3 P 9 - 4p r + -p q p T + p 2 T ) f + 

2 / 3 
- (6 + p q - 6 p T + 4 p q p T - 6 p 2 T + p q p\ + 6 pi) — + 

Pq 

3(l- Pr ) 2 (-l-2 P(? + 2p r -2p,p T -p2)Z! + 4(l-p r ) 4 4] + 

Pq Pq 

G b [2x+^x 2 + Ap q - ^p r + (-2 + 3p g - ^p 2 q -2p T + 5p qPr )f + 



I P 2 q + 4 Pr + 8 Pq pr ~ | p\ Pr - 2 pi - 10 p q p\ ) — 

O O p,j 



where 



and 



5 25 f 3 f 4 

•l + Pr)(3 + -p g -8p r + -p 9 p r + 5p^ + 5(l-p r ) 3 H}, (3) 

O O p g Pq 



1 + p T — X 



,2 G^mf . . 

^fel mo -3 - ( 5 ) 



192 7T 3 

The quantities and are the following hadronic matrix elements: 



1 <H b \bn 2 b\H b > /ml (6) 



2M Hb 
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Gb = ^ir~ < H »\ b \ ^G^\H b > /ml (7) 

where v is the 4- velocity of the decaying heavy meson = d^ — igA^, = A^T a is the 

gluon field in matrix representation, 7? is the residual part of the heavy quark momentum 
arising from its interaction with the light degrees of freedom and = gG a ^ u T a is the gluon 
strength tensor. 

The energy distribution @ spans in x from x = 2 y/pT- to x = 1 + p T — p q . It includes the 
nonperturbative corrections - terms proportional to Kb and Gb- The part of eq.(H) without 



corrections coincides with the electron spectrum in p decay with a massive from ref. |L3 



In the limit p q — > we encounter the familiar end-point singularities of the lepton spectrum. 
To get the limit right we make the following substitutions: 



r ^ 5(i + Pr -x) 

Pq n-l 

and 

fn 5 >( 1 + pr _ x ) 



p\ (n - 1) (n 



n > 1, (8) 



n > 2 (9) 



and only then take p q to zero. To see that this procedure gives the right limit, it is sufficient 
to compare the results of integration of the both sides of (H) and @ in the limits 2 y/pT- < 
x < 1 + p T — p q , multiplied by an arbitrary integrable function and check that the both sides 
are equal to each other. The distribution for p q —>■ takes the form: 

p9_>0 - • / x 2 -4p T {(6 + 6p T )x-4x 2 -8p T + 



r dx 

^ r 28p 1 4p r 2 4p3 p^ i 0x 2 
6 ^ ~1T~ + <5 (1 + p r - ^) (— - + — 2 p r + — — ) — ] + 

_ , 28p T e/1 w 11 _ 2 5p;_. , 10x 2 ,, 

[ ^ + 5(1 + Pr - x) {-— + 9 p r - 7p 2 + -p) + Ax + — ]}. (10) 



Now we can integrate distribution (|3|) to get the decay width including 1 / m 2 . corrections: 
F = l^l 2 tI^ 1 + ) z o( Pg ,Pr) -2G bZl ( Pq ,p T ) }, (11) 



where: 



Zo(p„Pt) = v^(l - 7p, - 7p* + p^ - 7p T - 7p= + p 3 T + p,p r (12 - 7p, - 7p r ) ) + 

i±^ + i2p; ( i-pj)iogii±^ 

1 - Vg) [l-V T 



>q < rq 

12 P l (1 - Pi) log + 12 P^ (1 - Pi) log ^4 (12) 



*1 (P 9 > Pr ) = ( (1 - Pq f + (1 - p T ) 3 ~ 1 - Pq Pt (4 - 7 p ? 

;i + 0(i + 



7Pr)) + 



12 plpl log 



«g)(l 



(13) 
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A = A(l, p q , p r ) = 1 + p\ + p 2 T - 2p q - 2p T - 2p q p T , 
v q and v T are the maximal velocities of the quark and the r-lepton produced in the decay: 

y/X y/X .... 

V r = 1 — — -■ ( 14 ) 



V, 



1 + Pq~ Pr 1 - Pq + Pr 

There exists a simple relation between the two functions z (p q ,p T ) and zi(p q , p T ) || : 
/ x dz (pq,p T ) dz (p q ,p T ) 

Zl{Pq,Pr) = ~2Pq -J^ 2p r — +AZo{pq,pr). (15) 

Width (|TT|) is symmetrical function of p q and p T and therefore its limit when p q — > 



r,^o = \V qb \ 2 S4{(1 + ) C 1 " 8 Pr + 8 Pi - Pr - 12 Pi log p. ) - 2 (1 - p T ) 4 } 



(16) 

is the same as limit when p T — > with substitution p q ^ p T . 



III. NUMERICAL ESTIMATES AND EXPERIMENTAL PREDICTIONS 

To make numerical estimates in this section we will consider the I?-meson decays B — > 
tz/JT. We choose the following values for the parameters entering expressions (|3|), fllCf ), ([TT]) 
and (|l6l) (see discussion in and fill): m b = 4.8 GeV, m c = 1.4 GeV, m u = 0, m r = 1.78 
GeV, Gb = 0.02 and K b = 0.017 (note that for A b we have = and the corrections are 
much smaller). 

For the width of the decay we see that the nonperturbative corrections are negative. 
For b — > mtz/ transitions they decrease the width by 6% of its perturbative value, while 
for the b — > ctv case they decrease it by 10% (note that for the massless lepton in the 
final state these numbers are 4% and 5% correspondingly). The quantity unsensitive to the 
uncertainties of known values of V qb is Y(b — > rvX)/T(b — > euX): 



r. 



Tjb^TlsX) = Z ( Pq ,p T ) ( Zl(pq,p T ) Zjjp^O) 

"~Y{b^evX) z ( Pq ,0) [ b[ z ( Pq ,p T ) z (p q} 0) n - 



The first factor of last equation describes the phase space suppression while the second 
contains the nonperturbative corrections. The corrections reduce r c by 4% and r u by 2% 
of their perturbative values. Note that in the quantity r c the relative contribution of the 
corrections is almost independent of the uncertainties in the value of c-quark mass. 

The obtained energy distributions @ and (|ToD could be applied to the decays involving 
the transitions b — » c and b — > u. The analysis of applicability of the distributions was made 



in P] and | 11[| . There it was shown that the proper quantity to confront with experiment is 



7 (x) = J x dx' jr^r, 2^p7 < x < 1 + p T - p q . (18) 

This quantity does not contain the end-point singularities and is suitable for direct com- 
parison with the experimental data for x not too close to its maximal value (so that the 
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operator product expansion is still valid and we are not in the resonance region). As we 
mentioned above, the hadronic kinematical region is different from the quark one. There- 
fore to compare our results with experiment the range of integration of the experimental 
distribution should include the window between the the quark and hadronic boundaries 
1 + p T — pg < x < Ms/mi, + p T — p q . Correspondingly, the reliable prediction can only 
be made for 2 v /p7 < x < x max , where x max = 1 + p T — p q — (Mb — rrib) jm^. For w-quark 
x max ~ 1.05, for c-quark x max ~ 0.95. 

The lepton energy spectrum is plotted on the Fig.l for b — > crv and on the Fig. 2 for 
b — > utv. The delta-functions of eq.(TD) are not shown on the graph. For comparison, on 
the same plots we show the energy distributions for electrons in b — ► cev and b —>■ uev 
transitions correspondingly (p T = 0). 

The function 7(x) is plotted on Fig. 3 for the case b — > ctv and on Fig. 4 for b — > utv. The 
solid line shows 7(x) with the nonperturbative corrections while the dashed line - without 
them. 

Although the functions are plotted for the whole range of x, 2^fp^ < x < 1 + p T — p q , we 
can only trust the graphs for x < x max . 



IV. CONCLUSIONS 

We investigated the nonperturbative corrections up to the order to semileptonic 

decays H b —* tvX and found their contribution to the r-lepton energy spectrum and the 
width of the decay. In the case of 5-meson decays the corrections could be up to 10% of 
the width. Unfortunately at present time the experimental measurements are not accurate 
enough in order to compare our results with the data. 
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APPENDIX: THE PHASE SPACE OF THE DECAY 

In this Appendix we briefly outline the derivation of the Lorentz Invariant Phase Space 
(LIPS) for the case of the inclusive semileptonic decay — > tvX . Let us introduce 
the invariant kinematical variables in the following way. Let P be the 4-momentum of the 
decaying particle, so that Mfj = P 2 ; q = p T + p p is the 4-momentum of the lepton pair, p T 
and pp are the momentums of the emitted charged lepton and neutrino correspondingly. We 
introduce m 2 x = p 2 x , - the invariant mass squared of the born hadronic state and m 2 Xv = p 2 X yi 
- the combined mass of the hadrons and v. We have the following relations: p x = P — q 



and p XD = P — p T . The introduced invariant variables are related to the ones used in |TT 
in the following way: 

m 2 x = M 2 +q 2 -2M Hb q , (Al) 
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m 2 Xp = M 2 +p 2 T -2M Hb E T . (A2) 



The inclusive decay rate in the normalization of [11] is given by the following expression: 



T = - J m2 XmaX dm\ \M(m 2 XP} m 2 x ,q 2 )\ 2 d<S>(P lPx ,p T} p p ) } (A3) 

X min 

where \A4(m Xp ,m x ,q 2 )\ 2 is the matrix element which describes the transition into the 
hadronic states with mass m x , it is given by the expressions (2)-(8) in The lower 

boundary of m 2 x is given by the mass squared of the lightest final hadronic state possible 
in the decay (D-meson mass for b — > c and pion mass for b — > u transitions ), while the 
upper boundary of m 2 x is (M^ — m T ) 2 . The d&(P,p x ,p T ,p p ) is the three particle LIPS for 
a particle with the 4- momentum P going into three particles with 4-momentums Px,Pt and 
p p . It is given by the formulae: 

dHP, Px ,p T ,p p ) = (2 vr) 4 5\P -px-pr- Pp ) 6(p x ) 9{p T ) 6(p P ) p* (A4) 

{lit) 6 (Air)- 3 [Ztc)^ 

where 6(a) denotes 6 (do), and is the zero component of the 4- vector a. For the three- 
particle phase space one can write: 

/d ijv^ - 
—^d^(P,p Xp ,p T )d^(p X p,Px,Pp)S(m XP -p 2 XP ) (A5) 
Z 7T 

where c?$ (a, b, c) is a two particle LIPS for a particle with 4-momentum a going into particles 
with 4-momentums b and c. The equation (]A5|) gives the decomposition of a 3-particle LIPS 
into 2-particle LIPSes. 

In order to introduce a new variable q 2 into the phase space integral one can multiply 
the equation (|A5|) by '1': 

J dq 2 5(q 2 -(p p + p T ) 2 )6(q ) = l. (A6) 

Performing the integrations over all variables except for m 2 x ,q 2 and m Xp , which amounts 
to integrating the delta functions in LIPSes along with determining the conditions for the 
delta functions to be non-zeros, we arrive to the following formulae: 

-(M Hb -m T ) 2 pm Xf> fq?na X ("i 2 x ,m 2 x -) 



512tt 2 M 



f\ al H b — m T) f m X- t'lraax\ ra X' X-> 

— / 2 dm\ P / dm 2 x dq 2 \M(m 2 XP ,m 2 x , q 2 )\ 2 , (A7) 

H b J K Jm d 2 J <L in {m x ,m xv ) 



where 



1 Ml -m 2 \A {Ml , ml, m\ P ) (m 2 P - m 2 ) 

2/2 2 \ 1 r i\*2 i 2 2,2 H b " 1 t 2 V v H ^ XuJ K Xu x, ~ 

<lmin( m Xi m Xu) = nl M H h + m X ~ m Xu + ™t ±o m X ~ 



Z m Xp m Xp 



(A8) 



2 x l r , ,2 2 2 2 M H h - m2 r 2 \^ K M k> m r> m Xu) {™Xv ~ ™l 

^x,m XP ) = -[M Hh +m x -m XP + m T ±-„ m x + 5 



(A9) 
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and 

X(a,b,c) = a 2 + b 2 + c 2 - 2 ab - 2 ac - 2bc. (A10) 

The physical meaning of the function X(a, b, c) lies in its relation to the square of the spatial 
momentum p* 2 of particles born in a two - body decay in the center of mass reference frame, 
for example: 



PXu-Pr ~ 4M 2 • l Aii J 



The fully differential distribution (|I|) in the invariant variables looks as follows: 
dT 1 



dm 2 x dm 2 XP dq 2 512 vr 2 M 2 H 



\M(m 2 x ,m 2 Xp ,dq 2 )\ 2 . (A12) 



S 
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FIGURES 



FIG. 1. The energy spectrum @ of r is plotted for b — > ctv transitions. The solid line shows 
the distribution with the nonperturbative corrections, while the dashed line - without them. For 
comparison, on the same plot we show the electron energy distribution for b — ► cev transitions 
(p T = 0). The graph can only be trusted for x < x max ~ 0.95. 

FIG. 2. The energy spectrum (|Tc|) of r is plotted for b — ► utu transitions (p q = 0). The solid 
line shows the distribution with the nonperturbative corrections, while the dashed line - without 
them. For comparison, on the same plot we show the electron energy distribution for b — ► ueu 
transitions (p T = 0). The graph can only be trusted for x < x max ~ 1.05. 

FIG. 3. The function 7(2:) plotted for the case b — ► ctv. The solid line shows 7(2;) with the 
nonperturbative corrections while the dashed line - without them. The graph can only be trusted 
for x < x max ~ 0.95. 

FIG. 4. The function j(x) for the case b — ► utu (p q = 0). The solid line shows 7(2;) with the 
nonperturbative corrections while the dashed line - without them. The graph can only be trusted 
for x < x max ~ 1.05. 
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